First-principles LDA+ U and GGA+ U study of neptunium dioxide 
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We have performed a systematic first-principles investigation to calculate the electronic structures, 
mechanical properties, and phonon dispersion curves of NpC>2. The local density approximation+[/ 
and the generalized gradient approximation+C/ formalisms have been used to account for the strong 
on-site Coulomb repulsion among the localized Np 5/ electrons. By choosing the Hubbard U pa- 
rameter around 4 eV, the orbital occupancy characters of Np 5/ and O 2p are in good agreement 
with recent experiments [J. Nucl. Mater. 389, 470 (2009)]. Comparing with our previous study of 
TI1O2, we note that stronger covalency exists in Np02 due to the more localization behavior of 5/ 
electrons of Np in line with the localization-delocalization trend exhibited by the actinides series. 
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I. INTRODUCTION 

Neptunium-based oxides have been extensively inves- 
tigated by experiments P, [H, H, 0, H, [1] and theoretical 
calculations 0, _§, [__, []J_] due to the technological impor- 
tance that neptunium is one of the important long-lived 
actinides which accumulate in high level waste during 
the conventional nuclear fuel cycle. Among neptunium 
oxides, Np02 has attracted much more attention because 
of its more stable thermodynamical properties and inter- 
esting subtle non-collinear magnetic structure [l], @, H, 0] • 
Np02 is reported to can maintain stable under high tem- 
perature up to 1400 K [1]. Seibert et al. systematically 
studied the reaction of neptunium with molecular and 
atomic oxygen as well as the the formation and stability 
of surface oxides; they reported that high neptunium ox- 
ides can only exist as surface phases not as stable bulk 
phases [6]. As for the electronic structures of actinidc 
dioxides (AO2), taking (A=U, Np or Pu)C"2 for exam- 
ple, the 5/ orbitals are partially occupied, accompanying 
the competition between the quantum process of local- 
ization and derealization of A 5/ electrons, which leads 
to many complex behaviors. This also makes the accu- 
rate description of electronic structure of AnC>2 series to 
be difficult achieved. 

Conventional density functional theory (DFT) that 
apply the local density approximation (LDA) or gen- 
eralized gradient approximation (GGA) underestimates 
the strong on-site Coulomb repulsion of the 5/ electrons 
and, consequently, describes actinide dioxides as incor- 
rect ferromagnetic (FM) metals instead of antiferromag- 
netic (AFM) Mott insulators. One more promising way 
to improve the drawback is LDA+U or GGA+U ap- 
proach, in which the underestimation of the strong in- 
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traatomic Coulomb interaction is corrected by the Hub- 
bard U parameter. Recently, the electronic structures, 
mechanical properties, and high- pressure behaviors of 
UO2 and PuC>2 have been correctly re pro duced using 
LDA+U or GGA+U calculations @, UMl M, Q El • 
The insulator character of NpC>2 is established by the 
experiments d, However, to our knowledge, a sys- 

tematical theoretical investigation of electronic struc- 
ture and mechanical properties for Np02 is still lack- 
ing. Consequently, based on the good performance of 
LDA/GGA+[/ approaches in describing the electronic 
structure of the strong-correlation systems, we carried 
out the present study of Np0 2 . 

In the present study, the lattice parameter, elec- 
tronic structure, mechanical features, and thermody- 
namic properties of Np02 are calculated by employing 
the LDA+f/ and GGA+U schemes due to Dudarev et 
al. We discuss how the choice of U as well as the 

choice of exchange-correlation potential, i.e., the LDA 
or the GGA, affect those properties. Our results show 
that the pure LDA or GGA fails to give the accurate 
lattice parameter and correct electronic structure, while 
the LDA+ U and GGA+ U schemes can effectively rem- 
edy these failures. The rest of this paper is arranged as 
follows. In Sec. II the computational method is briefly 
described. In Sec. Ill we present and discuss our results. 
In Sec. IV we summarize the conclusions of this work. 



II. COMPUTATIONAL METHODS 

Our total energy calculations are carried out by em- 
ploying the plane-wave basis pseudopotential method 
as implemented in Vienna ab initio simulation package 
(VASP) [Hj]. The exchange and correlation effects are 
described by the DFT within LDA and GGA [H, HO). 
The projected augmented wave (PAW) method of Blochl 
PH is implemented in VASP with the frozen-core approx- 
imation. Electron wave function is expanded in plane 
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waves up to a cutoff energy of 500 eV, and all atoms are 
fully relaxed until the Hellmann-Fcynman forces become 
less than 0.02 eV/A. A 9x9x9 Monkhorst-Pack [n| k 
point-mesh in the full wedge of the Brillouin zone is used 
for fluorite NpC>2. The neptunium 6s 2 7s 2 6p 6 6d 2 5/ 3 and 
oxygen 2s 2 2p 4 electrons are treated as valence electrons. 
The strong on-site Coulomb repulsion among the local- 
ized Np 5/ electrons is described by using the formalism 
formulated by Dudarev et al. [T3|- In this scheme, the 
total LDA (GGA) energy functional is of the form 



LDA(GGA)+U 



-^LDA(GGA) 
U - J 



[Trp CT -Tr(pV)], (1) 
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where p a is the density matrix of / states with spin 
cr, while U and J are the spherically averaged screened 
Coulomb energy and the exchange energy, respectively. 
In this work, the Coulomb U is treated as one variable, 
while the parameter J is set to 0.6 eV. Since only the 
difference between U and J is meaningful in Dudarev's 
approach, therefore, we label them as one single param- 
eter U for simplicity. 

Both spin-unpolarized and spin-polarized calculations 
are performed in this study. Compared with FM and 
AFM phases, the nonmagnetic (NM) phase is not en- 
ergetically favorable both in the LDA+ U and GGA+ U 
formalisms. Therefore, the results of NM are not pre- 
sented in the following. The dependence of the total 
energy (per formula unit at respective optimum geome- 
tries) on U for both FM and AFM phases within the 
LDA+ U and GGA+ U formalisms is shown in Fig. [T] 
On the whole, the total energy of FM phase is lower 
than that of AFM phase no matter in LDA+ U scheme 
or GGA+f/ scheme. However, as shown in Fig. 1, it is 
clearly that the difference can be negligible, especially as 
the increasing of the U parameter. The total-energy dif- 
ferences (Efm—Eafm) within the LDA+U* at U=4 and 
6 eV are —0.005 eV, which is close to previous calculation 
@ (-0.004 eV for [7=4.25 eV), and -0.002 eV, respec- 
tively. We stress that our approach does not include the 
spin-orbit coupling (SOC); a detailed explanation on this 
issue can be found in Ref. [1| . Both FM and AFM re- 
sults will be presented in the following analysis. 

In this study, the theoretical equilibrium volume, bulk 
modulus -B, and pressure derivative of the bulk mod- 
ulus B' are obtained by fitting the third-order Birch- 
Murnaghan equation of state (EOS) [HI- The elastic 
constants, various moduli, and Poisson's ratio v are cal- 
culated using the same method of our previous work [25| . 
Note that the bulk modulus B obtained by these two 
approaches are in good agreement, indicating that our 
calculations are self-consistent. 



FIG. 1: (Color online) Dependence of the total energies (per 
formula unit) on U for FM and AFM Np0 2 . 



III. RESULTS 
A. Atomic and electronic structures of Np0 2 

At room temperature and zero pressure conditions, the 
stoichiometric neptunium dioxide crystallize in a CaF2- 
like ionic structure with space group FmSm (No. 225), 
which is the stable structure of all actinide dioxides. Its 
cubic unit cell (defined by the lattice parameter ao) is 
composed of four Np02 formula units with the neptu- 
nium atoms and the oxygen atoms in 4a and 8c sites, re- 
spectively. At higher pressure, ^33 GPa experimentally, 
Benedict et al. [§} reported Np02 undergoes a phase 
transition to an orthorhombic Cmcm (No. 63) structure. 
In this study, we only focus our sight on the Fm3>m Np02. 

The present optimized lattice constant ao and bulk 
modulus B, obtained by fitting the Murnaghan equation 
of state for AFM Np02, are presented in Fig. O For 
comparison, the experimental values of ao (Ref. [2|) and 
B (Ref. [3J]) are also shown. Note that the experimental 
value of ao (5.42 A) Q is the fitted value at K rather 
than at ambient condition. In the overall view, the de- 
pendence of the lattice parameter ao of Np02 on U is 
similar to our previous study of Pu02 [IH- For the pure 
DFT calculations (Z7=0), as shown in Fig. Ufa), both 
LDA and GGA underestimate the lattice parameter with 
respect to the experimental value. This trend is more 
evident for LDA due to its overbinding character. After 
turning on the Hubbard U, for the LDA+ U approach, 
although the lattice parameter is still underestimated in 
a wide range of U, the calculated ao improves upon the 
pure LDA by steadily increasing its amplitude with U. 
Actually, at a typical value U =4 eV, the LDA+ U gives 
ao=5.40 A which is very close to the experimental value. 
On the other hand, the GGA+ U enlarges the underbuy- 
ing effect with increasing Hubbard U. As a comparison, 
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FIG. 3: (Color online) The total DOS for the Np0 2 AFM 
phase computed in the (a) LDA+ U and (b) GGA+ U for- 
malisms with four selective values of U. The projected DOSs 
for the Np 5/ and O 2p orbitals are also shown. The Fermi 
energy level is set at zero. 



FIG. 2: (Color online) Dependence of the (a) lattice parame- 
ter and (b) bulk modulus on U for AFM Np0 2 . 



at U=A eV, the GGA+ U gives a =5.50 A. Totally speak- 
ing, both the LDA+ U and GGA+ U results of the lattice 
parameter for the Np02 AFM phase are comparable with 
the experimental value at U around 4 eV. We have also 
calculated the equilibrium lattice parameter for the FM 
phase for Np02- The tendency of arj with U is similar to 
that of the AFM phase. Note that the screened Coulomb 
hybrid functional [7| and self-interaction corrected local 
spin-density (SIC-LSD) [10] calculations predicted the 
lattice parameter to be 5.42 and 5.46 A respectively. 

As for the dependence of bulk modulus B of AFM 
Np02 on U shown in Fig. Efb), it is clear that the 
LDA results (220-228 GPa) are always higher than the 
GGA results (188-200 GPa), which is due to the above 
mentioned overbinding effect of the LDA approach. At 
a typical value (7=4 eV, the LDA+U and GGA+U 
give 5=228 and 199 GPa, respectively. Apparently, the 
GGA+ U approach with increasing U gives more close 
values to the experimental data of 5=200 GP a [p . Note 
that the recent LDA+ U 3 and SIC-LSD [ltf calcu- 
lations predicted the bulk modulus to be 228 and 217 
GPa for AFM Np02, respectively. Overall, comparing 
with the experimental data and recent theoretical re- 
sults, the accuracy of our atomic-structure prediction for 
AFM NpC>2 is quite satisfactory by tuning the effective 
Hubbard parameter U in a range of 3-4 eV within the 
LDA/GGA+ U approaches, which supplies the safeguard 
for our following study of electronic structure and me- 
chanical properties of Np02 . 

Almost all the macroscopical properties of materials, 



such as hardness, elasticity, and conductivity, originate 
from their electronic structure properties as well as chem- 
ical bonding nature. Therefore, a correct description 
of electronic structure is necessary for further investi- 
gating these properties. In the following, we will show 
the dramatic improvement brought by the LDA/GGA+ U 
schemes compared to the pure ones as describing the elec- 
tronic structure properties of Np02- The total density of 
states (DOS) as well as the projected DOS for the Np 
5/ and O 2p orbitals at four selective values of U within 
LDA/GGA+f/ formalisms are shown in Fig. [3J With- 
out accounting for the on-site Coulomb repulsion ( U =0), 
one can see that both LDA and GGA methods predict 
an incorrect metallic ground state with nonzero occu- 
pation of Np 5/ states at Ep. After switching on U, 
as shown in Fig. [3J the Np 5/ bands begin to split at 
Ep and tend to open a Mott-type gap E g . The ampli- 
tude of this insulating gap increases with increasing U 
(see Fig. 4). Overall, the LDA+U and GGA+U give 
an equivalent description of the one-electron behaviors 
in a wide range of U. At a typical value U =4 eV, one 
can see from Fig. [3] that the occupied DOS near the 
Fermi level is featured by the three well-resolved peaks. 
The narrow one near —1.0 eV is principally Np 5/ in 
character, while the other two broad peaks respectively 
near —2.0 and —4.0 eV are mostly O 2p. These well- 
separated orbital peaks have been observed in the recent 
ultra-violet photoelectron spectroscopy (UPS) measure- 
ment Q . In addition, The width of the Np 5/ valence 
band is 1 eV, consistent with the hybrid functional cal- 
culation 0] and experimental data @| . The O 2p valence 
band width is of 4.5 eV, also in accord with the experi- 
mental data (4 eV) and previous hybrid functional result 
(4.5 eV). The Mott gap opened at [7=4 eV is of 2.2 (2.4) 
eV within the LDA(GGA)+ U scheme, see Table H Pre- 
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vious hybrid functional calculations result in a larger E g 
by ^0.8 eV [7]. The calculated amplitude of local spin 
moment is ~3.1^b (per Np atom) in both AFM and FM 
phases and within the two DFT+{7 schemes. On the 
whole, although the pure LDA and GGA fail to depict 
the electronic structure, especially the insulating nature 
and the occupied-state character of NpC>2, our present re- 
sults show that by tuning the effective Hubbard parame- 
ter in a reasonable range, the LDA/GGA+ U approaches 
can prominently improve upon the pure LDA/GGA cal- 
culations and, thus, can provide a satisfactory qualitative 
electronic structure description comparable with the ex- 
periments. By further increasing U to 6 eV, one can see 
that the peak near —1.0 eV becomes weak and is mostly 
O 2p. Meanwhile, there develops a prominent hybridiza- 
tion between Np 5/ and O 2p orbitals. This picture of 
DOS is no longer valid since the peak near the valence 
band maximum has been experimentally verified 0] to be 
due to the Np 5/ contribution. Thus, the overestimation 
of U will compel Np 5/ orbitals to be even more local- 
ized, resulting in their unphysical alignment with the O 
2p orbitals. 

In order to understand the chemical bonding nature 
of ground state Np02, we plot in Fig. [5] the valence 
charge density map of the (110) plane. Evidently, the 
charge density around Np and O ions are all near spher- 
ical distribution with slightly deformed toward the di- 
rection to their nearest neighboring atoms. Moreover, 
the charge density around Np and O ions is high while 
there are almost no remaining valence charges in the 
large octahedral-hole interstitial region. This suggests 
that remarkable ionization of neptunium and oxygen ions 
and significant insulating property exist for AFM Np02 . 
We have also plotted the line charge density distribu- 
tion along the nearest Np-0 bond (not shown) and find 
that the minimum value of charge density is 0.51 e/A 3 . 



The distance between this minimum point and Np ion 
is 1.29 A. Thus, the ionic radius of Np in Np02 can 
be defined as 7"n p =1.29 A. After subtracting this value 
from the Np-0 bond length (2.34 A), the oxygen ionic 
radius ro=l-05 A can be obtained. The corresponding 
valency of Np and O in Np02 can be estimated by cal- 
culating the valence charges within the spheres of ionic 
radii. After integration, we find 11.386 electrons around 
Np and 6.233 electrons around O. Therefore, the valency 
of Np and O can be represented as Np 3 614+ and O 233_ 
respectively, which demonstrates prominent ionicity of 
the Np-0 bond. However, the ionicity of Np-0 bond is 
relatively weaker when compared with the Th-0 bond 
in Th0 2 HI]. On the other hand, the Np-0 bond has 
stronger covalency than the Th-0 bond due to that the 
minimum value 0.51 e/A 3 of charge density along the Np- 
O bond in Np02 is prominently larger than that along 
the Th-0 bond (0.45 e/A 3 ) in Th0 2 [25[. This find- 
ing is consistent with the fact that the 5/ orbitals are 
more localized in the former in line with the localization- 
delocalization trend demonstrated by the actinides series. 
Through analyzing the orbital-resolved partial densities 
of states (PDOS) for the AFM Np0 2 , we find that the 
ionicity of the Np-0 bond can be attributed to the charge 
transfer from Np Qd and 5/ states to O 2p states, while 
the covalency of the Np-0 bond is associated with the hy- 
bridization of O 2s and Np 6p states located from —20.73 
to —13.07 eV. Overall, the Np-0 bond can be described 
as a mixture of covalent and ionic components. 



B. Mechanical properties of NpC>2 

Our calculated elastic constants, various moduli, pres- 
sure derivative of the bulk modulus B , and Poisson's 
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TABLE I: Calculated elastic constants, various moduli, pressure derivative of the bulk modulus B , Poisson's ratio V, spin 
moments (fJ. mag .), and insulating band gap (E g ) for FmZm NpC>2 at GPa. For comparison, experimental values from Ref. 
are also listed. 



Alagnetism 


Method 


(GPa) 


(GPa) 


(GPa) 


B 
(GPa) 


B 


G 
(GPa) 


E 
(GPa) 


(1 


f-Lmag. 
\H-B ) 


E 

(eV) 


AFM 




405.1 


128.8 


78.9 


221 


4.1 


98.9 


258.2 


0.305 


2.678 


0.0 




LDA+ U ( 17=4) 


399.5 


145.5 


72.9 


230 


4.5 


91.2 


241.6 


0.325 


3.074 


2.2 




GGA+ U (17=0) 


358.7 


105.2 


53.5 


190 


4.2 


76.2 


201.7 


0.323 


2.903 


0.0 




GGA+ U ((7=4) 


363.6 


118.8 


57.4 


200 


4.3 


78.1 


207.5 


0.327 


3.092 


2.4 
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LDA+(7 (17=0) 


415.0 


148.7 


79.3 


237 


4.5 


97.7 


257.8 


0.319 


3.054 


0.0 




LDA+ U ((7=4) 


403.6 


143.4 


73.7 


230 


4.2 


92.7 


245.2 


0.323 


3.090 


2.1 




GGA+U ((7=0) 


372.7 


116.8 


70.9 


202 


4.4 


90.0 


235.1 


0.306 


3.104 


0.0 




GGA+(7 ((7=4) 


365.8 


117.8 


57.5 


200 


4.5 


78.6 


208.6 


0.327 


3.103 


2.4 



ratio v for Frr&m Np0 2 at GPa are collected in Table 
UJ Here, the moduli and Poisson's ratio are deduced from 
the calculated elastic constants using the same scheme 
as that in our previous study on isostructural TI1O2 [H| . 
Each derived bulk modulus B turns out to be very close 
to that obtained by fitting Murnaghan equation of state. 
This indicates that our calculations are consistent and re- 
liable. At a typical value of U =4 eV, LDA+ U gives bulk 
modulus of 230 GPa for both AFM and FM Np0 2 , while 
GGA+ U gives 200 GPa. Obviously, GGA+ U approach 
gives more closer value than LDA+ U with respect to the 
experimental value. Overall, the bulk modulus of Np0 2 
is bigger than that of Th0 2 [HI]. This originates from 
relative higher valence electron density and shorter bond 
distances of Np0 2 than those of Th0 2 . 

In addition, the hardness of Np0 2 is inve stig ated by 
using the approach raised by Simunek et al. [26j . In the 
case of two atoms 1 and 2 forming one bond of strength 
S12 in a unit cell of volume f2, the expression for hardness 
has the form [26] 

H = (C/n)b 12 s 12 e-^, (2) 

where 

S12 = \J (eie 2 )/(nin 2 (ii2), e; = Z^/n (3) 

and 

h = (^-^f = 1 - [y/M/(ei + e 2 )f (4) 
ei + e 2 

are the strength and ionicity of the chemical bond, re- 
spectively, and d\2 is the interatomic distance; (7=1550 
and cr=4 are constants. The radius Ti is chosen to make 
sure that the sphere centered at atoms i in a crystal 
contains exactly the valence electronic charge Zi. For 
fluorite structure Np0 2 , &i 2 =32 counts the interatomic 
bonds between atoms Np (1) and O (2) in the unit cell, 
711=8 and n 2 =4 are coordination numbers of atom Np 
and O, respectively, ri = 1.715 (A) and r 2 =1.000 (A) are 



the atomic radii for Np and O atoms, respectively, Zi = lh 
and i? 2 =6 are valence charge for Np and O atoms, re- 
spectively, di2=2.34 (A) is the interatomic distance, and 
f2=157.42 (A 3 ) is the volume of unit cell. Using Eqs. 
(2)-(4), we obtain si 2 =0.0967 and / 2 =0.0347. The hard- 
ness of Np0 2 at its ground-state fluorite structure is thus 
given by ff=26.5 (GPa). This indicates that the fluorite 
Np0 2 is a hard material and approaches to a superhard 
material (hardness > 40 GPa). The high hardness of this 
crystal can be understood from the dense crystal struc- 
ture, which results in high valence electron density and 
short bond distances. We notice that the hardness of 
Np0 2 is higher than that of Th0 2 (22.4 GPa) [H|. This 
also can be attributed to the higher electron density and 
shorter bond distance of Np0 2 compared with Th0 2 . 



C. Phonon dispersion curve of NpC>2 

Before phonon dispersion calculation, we have calcu- 
lated the Born effective charges and dielectric constants 
of Np0 2 because of their critical importance to correct 
the LO-TO splitting. The Born effective charge Z* is 
a measurement of the change in electronic polarization 
due to the ionic displacements. The form of Z* for a 
particular atom depends on the atomic position symme- 
try. For fluorite Np0 2 , the effective charge tensor for 
both Np and O that occupy the 4 a and 8 c Wyckhoff po- 
sitions respectively are isotropic. After calculation, the 
Born effective charges of Np and O ions for AFM (FM) 
Np0 2 are found to be Z£ c =5.24 (5.26) and Z£=-2.62 
(—2.63), respectively, within LDA+ U formalism with the 
choice of 17=4.0 eV. In addition, the macroscopic static 
dielectric tensor is also isotropic and our computed value 
of dielectric constant is 5.58 for the AFM phase and 
5.64 for the FM phase. 

Employing the Hellmann-Feynman theorem and the 
direct method, we have calculated the phonon curves 
along some high-symmetry directions in the Brillouin 
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FIG. 6: (Color online) Calculated phonon dispersion curves 
(left panel) and corresponding PDOS ( right panel) for Np02- 
Both the AFM and FM results are calculated within LDA+ U 
formalism with U =4 eV. In right panel, the orange and green 
lines stand for AFM and FM results, respectively, while solid 
and dashed lines are used to distinguish the PDOS of Np and 
O. 

zone (BZ), together with the phonon density of states. 
For the phonon dispersion calculation, we use the 
2x2x2 fee supercell containing 96 atoms and the 4x4x4 
Monkhorst-Pack fc-point mesh for the BZ integration. In 
order to calculate the Hcllmann-Feynman forces, we dis- 
place four atoms (two Np and two O atoms) from their 
equilibrium positions and the amplitude of all the dis- 
placements is 0.03 A. The calculated phonon dispersion 
curves along the T— X— K— T— L— X— W— L directions 
are displayed in Fig. [6] for both FM and AFM phases 
ofNp0 2 . 

To our knowledge, no experimental phonon frequency 
results have been published for Np02- In the CaF 2 - 
structure primitive cell, there are only three atoms (one 
Np and two O atoms). Therefore, nine phonon modes ex- 
ist in the dispersion relations. One can see that the gap 
between the optic modes and the acoustic branches is not 
so evident. But the LO-TO splitting at Y point is very 
significant since the inclusion of polarization effects. Due 



to the fact that neptunium is heavier than oxygen atom, 
the vibration frequency of neptunium atom is lower than 
that of oxygen atom. Therefore, the phonon density of 
states for Np0 2 can be viewed as two parts. One is the 
part lower than 5.7 (5.0) THz for AFM (FM) Np0 2 where 
the main contribution comes from the neptunium sublat- 
tice, while the other part higher than 5.7 (5.0) THz are 
dominated by the dynamics of the light oxygen atoms. In 
addition, one can see from Fig. [6] that the difference in 
the phonon dispersion between the FM and AFM phases 
of Np0 2 is neglectable in most regions of the BZ, except 
for one TO branch, which is a little lower for the FM 
phase than for the AFM phase. 

IV. CONCLUSION 

In conclusion, we have studied the structural, elec- 
tronic, mechanical, and thermodynamic properties of 
Np0 2 within the LDA+ U and GGA+ U formalisms. The 
atomic structure, including lattice parameters and bulk 
modulus, and the one-electron behaviors of 5/ state have 
been systematically investigated as a function of the ef- 
fective on-site Coulomb repulsion parameter U. By choos- 
ing the Hubbard U parameter around 4 eV within the 
LDA/GGA+f/ approaches, most of our calculated re- 
sults are in good agreement with the experiments. Based 
on this, we have further investigated the electronic struc- 
tures of Np0 2 and the Np-0 bond nature. The insulating 
ground state can be well reproduced and the Np-0 bond 
is found to be of a mixture of both the ionic and cova- 
lent characters. As for the mechanical property, we find 
that GGA+ U approach gives more closer value of bulk 
modulus than LDA+ U with respect to the experimen- 
tal value and the hardness of Np0 2 is calculated to be 
26.5 GPa. For the thermodynamic analysis, the Born 
effective charges and dielectric constants as well as the 
phonon dispersion curves have also been presented. 
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